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A modified adiabatic method, previously tested in positron- hydrogen 
scattering, has been applied to the case of positron- helium scattering below 
the positronium threshold (17.8 ev). Since the ratio of the positron's velocity 
to that of the atomic electrons is about the same as in the hydrogen test easel, 
similar accuracy is expected, the main source of error being due to the 
approximate atomic wave function used. Phase shifts for L = 0,1,2 are presented, 
and total and momentum-transfer cross-sections are computed. The latter dis- 
agree completely with the only experiments performed to date. 
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INTRODUCTION 

In a recent paper 1 , the author applied the exact second-order polarization 
potential of Dalgarno and Iynn 2 to the scattering of low-energy positrons by 
atomic hydrogen. Several objections may be brought against this procedure. 

First, the polarization potential used is adiabatic, and does not account 
explicitly for velocity- dependent effects. A comparison of the adiabatic 
results 1 with the detailed variational s-wave phase shifts of Schwartz 3 indicates 
that the polarization potential is too attractive, and a semi-empirical modifi- 
cation was introduced to achieve essentially exact agreement. This modification 
consisted in suppressing most of the short-range (monopole) parts of the 
potential. Besides the good agreement thereby achieved, the justification for 
such a modification was based on the fact that an adiabatic method is most 
likely to be in error near the atomic nucleus, where complicated correlations 
can occur. The monopole suppression used was energy independent, but nevertheless 
excellent s-wave phase shifts resulted over the whole energy range up to the 
positronium threshold (6.8 ev). 

Second, the possibility of virtual positronium formation is not explicitly 
included. Since the positronium atom is a large, loosely bound system, it as 
not too surprising that s-wave scattering is not seriously affected by this 
omission, but higher impact parameter (higher angular momentum) scattering 
should include more virtual positronium. This is borne out by recent L = 1 
and L ~ 2 lower bound calculations 4 , which indicate that the polarization 
potential method underestimates the attraction in these partial waves. It 
is apparently not possible to retain the good s-wave agreement and also fit 
the p- and d- wave phase shifts with a simple, L- independent local potential. 
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The situation for positron- hydrogen scattering is clear: the modified 

adiabatic method will underestimate the total cross-sections as the positronium 
threshold is approached, and will also underestimate the angular dependence 
of the scattering. 

The purpose of the present paper is to apply the polarization potential 
method to the elastic scattering of positrons by helium atoms. This problem 
is presently too difficult for rigorous solution, but is much more accessible 
to experiments. This is due to the relatively high threshold for positronium 
formation in helium ( 17-8 ev) compared to that of hydrogen, and also due to 
the fact that atomic helium is stable. It seems probable that energy resolved 
positron beams will soon be available 5 , and the helium scattering experiment 
should prove feasible. At the moment, one experimental cross-section value 
exists 6 , although the analysis 7 is rather intricate, since it involves the 
diffusion of a swarm of positrons through gaseous helium under the influence 
of an electric field. 

It is, thus, of interest to examine the modified adiabatic method as applied 
to the helium case. To do this with any hope of success, one must first verify 
that the velocity range of interest is appropriate. For the case of hydrogen, 
the mean kinetic energy of the atomic electron is 13 . 6 eV and that of the 
positron is less than 6.8 eV. For Helium, each electron has an average kinetic 
energy of 39«5 eV, while that of the positron is less than 17-8 eV. One can 
see that the positron moves somewhat slower relative to the electron velocity 
for helium than for hydrogen. Thus the semi- classical criterion for an adiabatic 
approximation is at least as well satisfied. 
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2 

To apply the Dalgarno-Lynn potential one must have a hydrogenic ground 
state for the target atom. Using a shielded, uncorrelated hydrogenic wave 
function for helium introduces the largest uncontrolled error into the present 
work, and it remains to he seen how much error in the scattering this approxi- 
mation causes. The only parameters in the theory are the effective nuclear 
charge, the monopole suppression factor and one more parameter which can he 
used to bring the static polarizability into agreement with experiment. 
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The Modified Adiabatic Calculation 

The non- relativistic Hamiltonian for the system consisting of a fixed 
helium nucleus, two electrons and one positron can be written as follows: 


H = H e + H,i + H x + V 

= ' i'w * % 1 
Hl1 " ufs&i 

V - z £ w [ nFfii]. 


( 1 ) 


Here, rx and r2 are the position vectors of the two electrons and x is the 
position vector of the positron. Atomic units are used with energies in 
Rydbergs, lengths in units of the Bohr radius, a© , and momenta in units of 

-l 

• 

The modified adiabatic method proposes to use a scattering wave function 
of the form 


( 2 ) 
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In this equation, 0(r ,r_ ) is the ground state wave function of the helium 

^_L *V«£- 

atom and satisfies 

[M K l2 -E.] $- 0 , (5) 

where E 0 is the ground state energy. The positron scattering function X^x) 
has appropriate boundary conditions, while G and G 0 are the first-order adiabatic 
correlation function and its monopole part, respectively. They satisfy the 
equations: 

(V - <v>) 4 > (W 

and 

[g„ 

l °> K - (kb) 

with 

<gC) = in 4>*Q 4> 

and 

< g ) - <p„y = o 

As in I, both the potential V and the correlation function G may be understood 
to be expanded in Legendre polynomials in the angles between and r x and r 2 . 

The parameter a thus measures the amount of P Q (monopole) correlation retained 
in ¥ k , and C will be discussed later. 
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Assuming that Y, has the correct form to satisfy the Schrodinger equation, 
one writes 9 

[H-£ a - l<][l l5) 

as the effective Schrodinger equation for the positron. 

Using Eqs(4a,b), this becomes 

[ H* - k’ + < v > v c < VG > + ^~ c ) <vO] \ C*)*o. (6) 

As in the hydrogen case discussed in I, the various terms can be identified 
as follows: 

<v> is the first-order (static) potential Vi(x), < VG > is the second- 
order (polarization) potential V 2 (x), and < VG 0 > is the second-order monopole 
polarization potential V 2o (x) . To obtain the phase shifts in this approximation 
it is necessaiy to make a partial-wave decomposition of X^, and solve the 
resulting equations in an entirely straightforward manner. 

Of course, the true eigenfunction 0 describing the helium ground state 
is not known exactly, and some reasonably simple approximate wave functions 
must be substituted, in order to evaluate these potentials. (This represents 
the only extra inaccuracy in the helium case not appearing in the hydrogen 
case.) In evaluating V x (x), the two- term analytic approximation to the best 
Hartree function was chosen. According to Byatt le , it yields a first order 
potential 



( 7 ) 
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where 

A = 2.4915, 6 = 5 .4W. 


A very similar potential was used without polarization "by Malik 11 . To evaluate 
the second-order potentials, an even simpler form is necessary, or else the 
equations for G and Go become intractable. Choosing the uncorrelated, shielded 
hydrogenic function 


<f>(S,n)= 7f'p e 




( 8 ) 


one can easily show that the second-order potentials can be obtained from 
the corresponding potentials in the positron- hydrogen case by the following 
rule: 



Uelii/m 


2 V 2 (?*) M 




( 9 ) 


The second- order potentials for hydrogen are known 2 , and are exhibited in I. 
It remains to give reasonable values for the three parameters a, 3 and C. 

If it had been possible to use an exact wave function for the helium 
atom (as was the case for hydrogen), the only parameter in the theory would 
be a, the amount of monopole distortion retained in the scattering function. 
This parameter ought to be between zero and unity, and the results of I favor 
a *** o. In the present calculation, both extreme values have been tried. 

In I the long-range behavior of the polarization potential was shown to 
be 


- 




3_ 

ZV 


do) 
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which shows the correct dipole polarizability, 9 /2. In the case of helium, 
the corresponding limit gives a polarizability of The measured 12 

polarizability is 1.376, and to obtain this value requires C/g 4 = 0.1529* 

Now, if g is determined by a variational calculation, the minimum energy is 
obtained for g = 27/16, and then C = 1.240. Another approach argues that the 
energy minimum is very flat and that it is preferable to set C = 1, from 
which it follows that g = I.5992. Both of these methods were examined 13 . 
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Results and Discussion 

Figs. 1, 2, and 3 show the energy dependence of the phase shift for 
L = 0, 1, 2 up to the threshold for positronium formation (E = 1.31) • The 
four curves represent the extreme values of the monopole suppression, and the 
two methods of adjusting the dipole polarizability discussed above. As the 
angular momentum increases these curves coalesce, and for L k J an analytic 
expression 14 involving only the dipole polarizability P is sufficiently accurate, 
i.e . 


* 7Y Pk S /(iL->){2L*i) (2L+$) 


(u) 


The total cross section 


o t. = 4 a (v-* 1 ) s.v s t te) 

T k l 

is shown in Fig. 4, for the same four cases, in units 

p 

of Too. In Fig. 5 the momentum- transfer cross-section is similarly displayed. 
This cross-section, defined as 


<r MT = 'fa 


(13) 


can be written (in units of Traf ) as 

cr' - 4 t S (iM 

°mt ~ k* X v J 


and is the . particular cross-section which is effective in determining the 
mean free path for swarm or diffusion measurements. Such a measurement has 
been madef and analyzed 7 under the assumptions that both cr^ and the 
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annihilation rate of free positions in helium are independent of energy, and 
that positronium formation is very rapid for any positrons attaining the 
threshold energy. The numerical value derived in this way is = (0.023 ± 
.006)rra^, and is shown in Fig. 5 "by the shaded band. It is seen that there 
is no agreement whatsoever between the experiment and any of the theoretical 
curves. We have repeated the diffusion analysis of Ref. 'I, allowing the cross- 
section to take the forms given by the present theory, but the disagreement 
remains. Note that we expect the L = 1 and L = 2 phase shifts to be under - 
estimated by the present calculation; for the diffusion analysis the higher 
energy range is most important. To estimate the effect on the momentum-transfer 
cross-section of increasing these two phase shifts one can differentiate Eq. (14) 
and replace the sine function by its argument. Then 


A 0”mt — 


(- Sc +3 5, - 2 : %) AS, f (-2S, ♦ 5S,-3§ s ) AS,] . 


(15) 


Inserting the values of the phase shifts at threshold, and assuming AS /6 « 0.5 
(which ip justified by comparison with the results of I) one finds that 
increasing 61 and 62 makes the agreement with experiment 6 * 7 even worse than 
before. The experiment and its analysis are not completely straightforward, 
and further experiments would seem to be desirable. 

One other calculation has been recently carried out 15 which deserves 
mention. It also employs an effective potential which is actually 

evaluated in the adiabatic approximation but, in common with similar adiabatic 
calculations in hydrogen, only a limited number of multipoles ( 0 , 1 , 2 ) were 
retained, instead of the complete polarization potential being used. Nevertheless, 
strong disagreement with experiment still occurs. 
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Figure Captions 

. Phase shifts in radians for s-wave scattering. The four sets of 
parameters discussed in the text are shown. The abscissa is in 
atomic units of momentum. 

• Phase shifts in radians for p-wave scattering. The four sets of 
parameters discussed in the text are shown. The abscissa here is 
in energy units. 

. Phase shifts in radians for d-wave scattering. The four different 
curves are very close, and are shown as a band. The abscissa again 
is in energy units, and the linear relation of Eq. (11) is shown to 
agree well. 

. Total cross-sections in units of tt^, including L = 0,1,2, 3 partial 
waves, for the four cases discussed. 

. Momentum-transfer cross-sections in units of na^, including L = 0,1, 
2,3 partial waves, for the four cases discussed. The "experimental 
result (Ref. 7) is shown by the shaded band. 
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